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PRELIMINARY MEASUREMENTS OF TURBULENCE AND TEMPERATURE 
FLUCTUATIONS BEHIND A HEATED GRID 


By A. L. Kistler, V. O'Brien, and S. Corrsin 
SUMMARY 


Preliminary measurements have been made of velocity and temperature 
fluctuations in the flow behind a heated grid in a uniform airstream. 
Temperature correlation shows a reasonable degree of isotropy, and the 
temperature fluctuations die out at large distances more slowly than the 
turbulence, as has been predicted theoretically under some strongly sin- 
plifying postulates. 


INTRODUCTION 


A comparison between the correlation equations for concomitant 
(incompressible isotropic velocity and temperature fields showed marked 
differences attributable to the fact that velocity is a vector while 


temperature is a scalar (refs. 1, 2, and پ١‎ .1 The relative decay rates 
were calculated under strongly simplifying assumptions, and the rela- 
tive "microscales” and relative wave numbers characterizing the very 
fine structures were deduced in terms of the Prandtl mumber of the fluid. 
Since the theoretical predictions are essentially conjectural, it is 
necessary to get at some facts through measurement. 


A hot grid with thermal mesh equal to momentum mesh is as simple 
as any, so this configuration was selected, with the expectation that 
approximate equality of integral scales would be obtained. 


The analysis of reference 1 suggests that the flow with equal inte- 
gral scales for velocity and temperature fields may not be the simplest 
case analytically, but could be the simplest for experimental realization. 


tas pointed out in references 1 and 2, the work applies equally well 
to isotropic turbulent mixing between two different components provided 
that the molecular mass transfer coefficient is nearly constant over the 
concentration range encountered. 
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This report describes the preliminary data. It is hoped that with 
improved grid and measuring procedure, less scattered results will be 
gotten in the measurements now in progress. 


The investigation was conducted at the Johns Hopkins University 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. 
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SYMBOLS 


one-dimensional power spectrum of u-fluctuations 


one-dimensional power spectrum of %$-fluctuations 


"longitudinal" double velocity correlation, 


u(x,y,z)u(x + £,y,2) 


u'(x)u'(x +£) 


"lateral" double velocity correlation, Sas eed + LE) 
UT 


wave number in x-direction 


spectral "cut off" for E- fk and E , respectively 
LUCE 8 (L 


integral scales of g, f, and m, respectively 

grid mesh size for momentum and heat (M = My = 1 in.) 
temperature correlation 

space interval 

time 

turbulent velocity fluctuation in x-direction 

Cartesian space coordinates; x is alineà with mean flow 
thermal diffusivity coefficient, K/ Cp, where K is 


thermal conductivity coefficient, p is density, and 
Cp is specific heat at constant pressure 
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1 space interval in y-direction, SAY 

3 temperature fluctuation 

As AS dissipation scales for velocity and temperature fluctua- 
tions, respectively 

۷ kinematic viscosity coefficient 

3 space interval in x-direction, =Ax 

o Prandti number, v/ Y 


EQUIPMENT AND PROCEDURES 


Approximately isotropic turbulence is produced at some distance 
behind a square-mesh biplane of grid of la inch round wooden dowels 
spaced 1 inch on centers. The grid is set at the upstream end of the 
working section of the wind tunnel sketched in figure l. 


As a first method of generating isotropic temperature fluctuations 
in this turbulence, a grid of Nichrome heating wires was introduced about 


1$ inches upstream of the dowel grid, with wires geometrically in phase, 


that is, alined so that their wakes passed over the rods (fig. 2). In 
this way the hot air passed through the turbulence-generating shear zones 
of the main grid, so that temperature and velocity fluctuations were pro- 
duced in the same local neighborhoods. 


Turbulence and temperature fluctuation measurements were made with 
the hot-wire anemometry equipment described in reference ¥. The experi- 
mental procedures for simultaneous velocity and temperature fluctuations 
are given in reference 5. ۱ 


The measurements were all made at a mean speed of 14 feet per sec- 
ond and a mean temperature rise (across the grid) of about 50 C. The 
actual heating wire temperatures were about 5009 C. The resulting tem- 
perature fluctuations were inconveniently small, but this insured neg- 
ligible influence of density variations upon the fully developed turbu- 
lent velocity field - as attested by both turbulence level and correlation 
measurements. 


Because of the extremely small temperature fluctuations it was nec- 
essary to use unusually high resistance hot-wires, 0.000l-inch platinum 
of about 90 chms resistence. This wire was set in the form of a V (with 
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plane perpendicular to mean flow) in order to reduce the error due to 
finite wire length. The high resistance necessitated correction for 
heating current fluctuations when the system was operated as an anemometer. 


As has been found before (ref. 6), the velocity field was not quite 
isotropic within the precision of measurement: v' ~0.9u' at x/M = 100; 
but this is apparently insufficient anisotropy to cause difficulties. A 
reasonable degree of isotropy in the thermal field is indicated by the 
approximate equality of سسجت‎ and lateral ieee correlation 
functions (fig. 5}. nn 


For convenience the spectrum was measured without heating, after it 
had been determined that the u'-levels and correlations were not appre- 


clably changed by the heating. 


The temperature spectrum was measured with "mixed" sensitivity, that 
is, wire temperature set for the same order of response to velocity and 
temperature. This was done because under resistance-thermometer condi- 
tions (heating current <1 milliampere) the signal was of the order of 
the amplifier noise level over most of the frequency range. 


With mixed sensitivity the hot-wire responds to both velocity and 
temperature spectra, which must be "separated." Simple superposition of 
the two energy spectra follows from the assumption that all "harmonics" 
of the two spectra are uncorrelated. Since uy = 0 in this field, such 
an assumption seems reasonable. Then the temperature spectrum is obtained 
by subtracting the (unheated grid) velocity spectrum from the mixed 
spectrum. 


A Hewlett-Packard model 300 A wave analyzer was used as variable 
frequency band pass filter. 


No correction has been made for wire length. 
EXPERIMENTAL RESULTS 


The velocity and temperature decay curves are plotted in figures 3 
and 4, respectively. A number of runs are given for each to indicate 
the repeatability of. the data as well as the equality of hot- and cold- 
grid flows (fig. 3) and the consistency of two difference measuring 
methods (fig. 4). The u'-values are a bit lower than those obtained 
behind a similar grid in a different wind tunnel (ref. 6). 


Figure 5 shows longitudinal and lateral temperature correlations 
at x/M = 69. Of course it was necessary to keep the downstream probe 
out of the wake of the upstream one, so the Ax coordinate is actually 
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inclined to the x-axis by about 5°. The two sets of points match fairly 
well, indicating that the temperature field is not far from isotropy. 


Figures 6 and 7 compare these curves with the corresponding velocity 
correlations. The close agreement of the longitudinal cases is especially 
noteworthy. 


The experimental scatter is too great to permit a determination of 
the dissipative scales ("microscales") by graphical fitting of parabolas 
to the vertices of these correlation data. However, the integral scales 
have been calculated: 


r= f g(n) an 


EL 


3 f(e) de (1) 
0 


A m(r) dr 


where f(r) and g(r) are the Karmdn-Howarth notation so that 


Lg 


f(E) = u(x,y,z)u(x + E,ÿY,2Z) 
u'(x)u'(x + رع‎ 


g(n) جرعاد(ترتگللاے‎ + az) 
(x) 


Also 


S(x,y, z) lx + &,y,2) 


(E) = ET + E) 
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and 


8(x,y,2)3(x,y + 152) 
8۶ (x) 


m(n) = 


are the longitudinal and lateral temperature correlation coefficient 
functions. 


The results are 


L = 0.40 in. 


Ly = 0.70 in. (2) 


ty = 0-66 in., laterally 
0.68 in., longitudinally 


In passing one notes that = = 0.57 £ 0.50, the value corresponding to 
E 


exactly isotropic turbulence. Since M, the momentum grid mesh, and 
Ms, the heat grid mesh, are both 1.0 inch, the values in equations (2) 
are also L/M, L1/M, and Lg/ Mg و‎ respectively. 


The time spectra of ð and u at x/M = 23 were determined as 
described in the previous section and are given in figure 8. These may 
be considered good approximations to the one-dimensional longitudinal 
space spectra, except perhaps for the lowest wave number range (refs. 7 
and 8). The sketched curves indicate a possible trend, but should not 
be taken seriously; within the experimental precision there is no sig- 
nificant difference. The circuit noise has been subtracted out for both 
sets of data. 


These spectra have been used for approximate calculation of the 
dissipative scales of velocity and temperature fields: 


ky “By (ky) akı (3)‏ ل“ 
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=- z ky “Ep (k1) akı (4) 


The results of graphical computation are 


in.‏ 0.26 = و۸ 


(5) 
À = 0.19 in. 


at x/M = 23. 


With the assumption that velocity and temperature fields are suf- 
ficiently near to isotropy A(x) and Ag(x) can be calculated from 


the intensity traverses, as indicated in the next section. The resulting 
faired curves are plotted in figure 9, with the above values spotted in. 


COMPARISON WITH THEORY 


Of course the theories of isotropic turbulence and of isotropic 
scalar fluctuations in isotropic turbulence are still unsolved; in fact 
they have not even been formulated in a determinate way. Nevertheless, 

a number of consequences of isotropy have been inferred, including some 
particular (and relatively simple) forms of the averaged differential 
equations. Of especial interest here are the so-called “decay equations." 
For isotropic fields decaying in time, there are (refs. 1 and 9) 


2 2 
du” س‎ — UT 
= 10v 2 (6) 
and 
2 2 
a 2 و‎ (7) 
Ag 


The actual experiments are done in steady state behind a grid with 
the various averaged quantities depending upon x instead of t. Taylor 
(ref. 9) has pointed out, in effect, that if these vary sufficiently 
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slowly with x equation (6) can be applied to the grid flow with the 
substitution x/U = t. The same holds for equation (7), so that 


du? 0 E 
os و‎ t (8) 
42 2 ۱ 
لے کے سے‎ nt (9) 
U Ag 


These two equations permit calculation of A(x) and Ag(x) from 


the experimental data on u'(x) and %'(x), as mentioned at the end of 
the previous section. The faired results are drawn in figure 9. The 
appreciable disagreement at x/M = 23 between A's calculated from 
decay and from spectra is indicative of inaccurate spectra (the less 
accurate of the two kinds of data) and/or of appreciable deviation from 
isotropy this close to the grid. 


The ratio A g/A varies almost linearly from 1.23 at x/M = 25 


to 1.58 at x/M= 100 (fig. 9). The only other measurement of this 
ratio was on the axis of a heated round turbulent jet (ref. 4), where 
the value was about 1.34. 


In reference 1 it is shown that for two special "kinds" of isotropic 
fields this ratio has the same simple value: 


)10( ك 


where o is the Prandtl number v/y. The two kinds of isotropic fields 
are > 


(a) Reynolds and Peclet numbers are so small that the convective 
effects are negligible for both heat and momentum. Therefore this is 
not really turbulence in the ordinary sense. 


(b) Reynolds and Peclet numbers are both very large and each field 
is assumed to have a completely "self-preserving" correlation function 
during decay. In this model, due to Von Kármán (ref. 10), the dissipation 


Tn fact, this requires perhaps more justification than given by 
Taylor or by Kármán and Howarth. This has been discussed at a seminar 
lecture at M.I.T., May 1952, and printed copies are available from the 
Aeronautics Department of The Johns Hopkins University: "Remarks on the 
Comparison of Grid-Produced Turbulence With the Theory of Isotropic 
Turbulence," by S. Corrsin. 


2N 
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and conduction terms of the two correlation equations are neglected, but 
decay is indirectly included by substitution of the decay equations into 
terms of the correlation equation. 


Equation (10) can also be deduced by an approach somewhat different 
from the above: 


Assume the Reynolds and Peclet numbers so high that the "inertial 
subranges " (vith رظ‎ (1) and Eg (ki) both ~ رت‎ 5) are sufficiently 
extensive to permit the approximations (ref. 2) 


O elsewhere 
(11) 
ai ka 219 
Eg(k) ~ ky 5/5 when gk, <k] < 4X 


O elsewhere 


The "cut off" wave number k, is just the inverse of Kolmogoroff's 
microscale, 


ke = کے‎ (12) 
and „gk, is the analogous thermal quantity (ref. 2) 


In reference 2 it was shown that 


Cesa” (14) 
ka 
Substituting equations (11), (13), and (14) into (3) and (4) and 


neglecting «2/3 relative to «2/3 and gk 2/3 relative to ske2/5, 
there results 
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۔ 2/3 2 
(=e) (15)‏ سے 


Tf one now restricts to fields in which the thermal large structure is 
roughly equal to the velocity large structure (as in these experiments), 
ko = Ko, and one arrives again at equation (10). 


In fact, none of the three sets of assumptions leading to equa- 
tion (10) is directly applicable to these experiments. The Reynolds | 
and Peclet numbers are not low enough for case (a) and not high enough 
for the other two cases. Furthermore, it is known that the correlation 
functions of velocity do not remain completely similar during decay 
(ref. 11). The important conjecture to be made is that a result encoun- 
tered under such divergent sets of assumptions may have approximate 
validity in a wide range of situations - perhaps including the actual 
one. 


Taking o = 0.72 for air, equation (10) gives 


DD 
>< 1.67 


which is in reasonable agreement with the downstream values computed 
from the decay data (fig. 9}. 


The Ag/A curve is just a particular way of displaying the decay 


data. Two other ways are instructive: Figure 10 shows the variation 
in Reynolds and Peclet numbers 


(16) 


The apparent constancy of Pag is evident. Figure 11 shows the relative 


decay rates of temperature fluctuations and velocity fluctuations. For 

small values of x/M they die out at about the same rate, but for larger 

values of x/M ð' dies out at a decidedly lower rate than does u'. 

For fully isotropic fields the ratio of equation (7) to equation (6) gives 5 
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NERF (17) 
dut/ut 5a \Ag 
and if equation (10) holds 
Eg- as) 


which is very close to the measured results at large values of x/M. 
This is just another way of showing that A j/A does approach the value 


in equation (10). 


In the absence of more complete data, no attempt has been made to 
compare the individual decay rates of 3' and wu’ with the theoretical 
predictions that can be obtained by a variety of postulates. 


CONCLUSIONS 


The principal conclusions to be drawn fram the foregoing results on 
measurements of turbulence and temperature fluctuations behind a heated 
grid are that: 


l. It is experimentally feasible to set up and measure concomitant 
velocity and temperature fluctuation fields downstream of a heated grid. 


2. The two fields are moderately isotropic sufficiently far 
downstream. 


3. Also, sufficiently far downstream, the dissipation scale ratio 
is approximately given by 
> کہ‎ 
À ر‎ 


a value deduced theoretically under widely divergent postulates, nore of 
which corresponds closely to the present experimental conditions . 
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L. As a consequence of conclusions 2 and 3, far enough from the 
grid the temperature fluctuations die out more slowly than velocity 
fluctuations: 


sx) ~ [a (x)] 3/5 


corresponding to equation (18). 


The Johns Hopkins University, 
Baltimore 18, Md., January 13, 1954. 
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heating mesh of double wires can be seen through wooden dowel grid. 
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Figure 3.- Turbulence decay. 
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Figure 4.~- Temperature fluctuation decay. 
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Figure 5,- Longitudinal and transversal temperature correlations. 
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Figure 6.- Transversal velocity and temperature correlations. 
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Figure 7.- Longitudinal velocity and temperature correlations. 
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Figure 8.- u and § spectra at x/M = 23.1. 
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Figure 9.- A's calculated from u' and 3' 
measured spectra). 
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Figure 11.- Relative decay rates of u' 
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